Hysteresis in the de Haas-van Alphen Effect 
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A hysteresis loop is observed for the first time in the de Ifaas-van Alphen (dHvA) effect of 
beryllium at low temperatures and quantizing magnetic field applied parallel to the hexagonal axis 
of the single crystal. The irreversible behavior of the magnetization occurs at the paramagnetic 
part of the dHvA period in conditions of Condon domain formation arising by strong enough dHvA 
amplitude. The resulting extremely nonlinear response to a very small modulation field offers the 
possibility to find in a simple way the Condon domain phase diagram. From a harmonic analysis, 
the shape and size of the hysteresis loop is constructed. 

PACS numbers: 75.45.+j, 71.70.Di, 75.60.-d 
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It is well known that the irreversibility of the magne- 
tization process by domain wall motion is due to energy 
barriers arising from a variety of defects inside magnetic 
material 1] . The resulting hysteresis effects in usual mag- 
netic substances, i.e. in substances where the atomic 
magnetic moments, due largely to the electron spin, are 
the reason of magnetism, have been investigated in detail 
in the past. However, in an applied magnetic field, the 
orbital motion of free electrons in metals leads also to 
magnetization. For Landau quantization 01 of the elec- 
tronic system, the oscillating magnetization known as the 
de Haas- van Alphen (dHvA) effect has been extensively 
studied in single crystals at low temperatures There 
is no experimental data reporting hysteresis in the dHvA 
effect apart from the trivial case caused by eddy currents 
when the applied magnetic field is varied fast enough 01 ■ 

If the dHvA amplitude becomes very large and com- 
parable to the oscillation period, i.e. ^qOM/OB > 1 
for magnetization M and magnetic induction B, a ther- 
modynamic instability arises at the paramagnetic (x — 
OM/dH > for susceptibility x and applied magnetic 
field H) part of every dHvA period, according to the 
Pippard-Shoenberg concept of high magnetic interac- 
tion _5] . For long samples with the applied field par- 
allel to the long axis (demagnetization factor n « 0) the 
instability is avoided by a jump of the magnetic induc- 
tion B between the stable states Bi and B2 at a criti- 
cal field He- For plate-like samples perpendicular to H 
(n « 1) i? = ^ioH is required so that the induction can 
not jump. In this case Condon domains arise with neigh- 
boring regions of respective inductions and B2 in the 
applied field interval Bi < /io-ff < ^2 0- For samples 
with < n < 1 this interval decreases proportionally to 
n. For samples of arbitrary shape there is an inhomoge- 
neous state - Condon domain state (CDS) with the same 
phases Bi and B2 ■ 



After the discovery of Condon domains [3,13 hysteresis 
in the dHvA effect needs to be considered again. Indeed, 



the CDS consists of two phases of different induction val- 
ues with a magnetization current in the domain walls. 
This needs usually extra energy. It was shown that the 
transition from the homogeneous state to the CDS is of 
first order At this phase transition one could expect, 
in principle, all phenomena like irreversibility, supercool- 
ing and hysteresis that exist at first order phase transi- 
tions, e.g. the liquid - gas transition. Naturally, Condon 
discussed these problems in his first paper on domains 
concluding that neither supercooling nor hysteresis had 
been observed in all reported data. Since then, the above 
phenomena were discussed in several papers [3, Il3 • De- 
spite Condon domains themselves have been observed 
and investigated experimentally [iTl IT^ IT^ IT^ . there 
is nevertheless up to now no experimental data proving 
hysteresis. 

In this Letter we present the first experimental ob- 
servation of hysteresis in beryllium single crystals. The 
hysteresis itself turns out to be very small. Therefore, 
several methods were used to prove the existence of hys- 
teresis. First the hysteresis effect is measured directly by 
Hall probes in DC fields, then a standard AC method 
is used with various modulation levels, frequencies, and 
magnetic field ramp rates. Finally, the DC hysteresis 
loop is reconstructed by assembling several higher har- 
monics of the AC pickup voltage. Moreover, it is shown 
that the Condon domain phase diagram can be measured 
directly from the response to sufficiently small AC fields. 
Finally, Plummer's [l5j strange and up to now only in- 
completely understood data is explained. 

Beryllium is to our knowledge the best metal to inves- 
tigate hysteresis effects related to Condon domain forma- 
tion. First, due to its cigar- like Fermi surface (small cur- 
vature at the maximal cross sections) the dHvA ampli- 
tude is very high for H || [0001]. Secondly, in this configu- 
ration two rather close dHvA frequencies 970 T and 940 T 
coexist leading to a beat in M and %. Because of this 
beat there is the unique possibility to change the dHvA 
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amplitude and the critical parameter % by a factor three 
at constant T by varying the magnetic field only very 
little. Thus, experimental conditions can be adjusted so 
that the transition to the inhomogeneous Condon do- 
main state occurs in a part of each beat period. The 
sample was cut from the same piece that was used ear- 
lier for the preparation of the plate-like sample in which 
Condon domains were first observed by /iSR The 
results shown here were measured on a rod-like sample 
of size 8x2x1 mm'^ with the long axis parallel to [0001]. 
The Dingle temperature was = 2.0 K. The measure- 
ments were made either directly by a micro Hall probe 
placed close to one end of the sample or by a compen- 
sated pickup coil using low frequencies of about 21 Hz 
and small modulation level (< 6 G). The experiments 
were carried out in a 10 T superconducting coil with ho- 
mogeneity better than 10 ppm in 1 cm"^. 

Fig. n shows Hall probe traces for an up and down 
sweep of the applied magnetic field around a dHvA antin- 
ode at 3.6 T and at T = 1.3 K. The hysteresis loop is 
very small and only visible by a zoom. A periodically 
arising induction difference 5B of about 3 G between the 
up and down sweep is measured at the steeper part, i.e. 
the paramagnetic part, of each dHvA period. The signal 
is about ten times higher than the noise level. The ap- 
plied field is measured here by another Hall probe, placed 
sufficiently far from the sample, as the superconducting 
solenoid has its own small hysteresis when the current is 
swept. The hysteresis effect is clearly observed in fig. ^ 
using DC Hall probes but rather under the most favor- 
able conditions at 1.3 K around a maximum of a mag- 
netization beat. The DC Hall technique is not sensitive 
enough to study hysteresis as a function of temperature 
and field. 

A standard AC modulation method with a compen- 
sated coil system used at very low modulation level h 
is much more sensitive to detect nonreversible magne- 
tization or hysteresis. This method is used to deter- 
mine exactly the point of its appearance. If the mod- 
ulation amplitude h is much smaller than the oscillation 
period Aff , the measured response corresponds in good 
approximation to the derivative, i.e. the susceptibility 
x{H) — dM/dH. If the modulation is further decreased, 
the result should not change, the measurements should 
only be more precise. 

Fig-Elshows a different behavior as a function of mod- 
ulation field amplitude h. All x{H) curves of Fig. \^&) 
and (b) are measured at low temperature where in the re- 
gions of each antinode hysteresis exists like it was shown 
on Fig.^ Fig.|2Ia) shows the normalized pickup voltage, 
i.e. response divided by the modulation level, measured 
by the AC method with 6 and 0.5 G modulation ampli- 
tude. In both cases h <^ AH so that one could expect 
identical curves. In fact, both curves are completely con- 
gruent except at the regions near the antinodes. Here, 
only the high modulation (6 G) level gives the expected 
result, which is the well known dHvA oscillation beat of 
beryllium. For comparison fig. |2Ib) shows x(-ff), calcu- 
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FIG. 1: Hysteresis loop observed by Hall probes in the para- 
magnetic part of the dHvA oscillations of beryllium. B{H) 
traces for an up and down sweep of the applied magnetic field 
around a dHvA antinode at 3.6 T and T — 1.3 K (scale on 
the left). The hysteresis loop is visible in the insert. The 
induction difference SB between these curves shows the value 
of the hysteresis (scale on the right). 



lated from B{H) curves measured with Hall probes like 
in fig. n without field modulation. For small modula- 
tion amplitude deep "notches" in the dHvA oscillation 
envelope are observed. The notches occur at magnetic 
fields where the dHvA amplitude is big enough that the 
above described instability arises. Instead of a further in- 
crease, after having crossed the critical point, the dHvA 
amplitude decreases at the paramagnetic (x > 0) part of 
the dHvA period only. The diamagnetic part does not 
change. One would expect that both M{H) and x{H) 
oscillate around zero. This shows that the small modula- 
tion measurements can not correspond to the real x{H) 
wave form. The amplitude decrease is absent at temper- 
atures T > 3.4 K and at the nodes of the dHvA beat 
where the single crystal is in the homogeneous state. 

The decrease of the dHvA amplitude at low modula- 
tion level is schematically shown on fig. Ola). If the mod- 
ulation is much bigger than the width of the hysteresis 
SH the AC response corresponds in good agreement to 
the slope, i.e. x{H)- As soon as the modulation becomes 
comparable to the hysteresis loop size the response de- 
creases . The effect was measured for several modula- 
tion levels at different temperatures and at different mag- 
netic fields with respect to the dHvA beat phase. The 
results are shown on fig. Olb) and (c). The modulation 
level where the dHvA amplitude decreases gives the hys- 
teresis width 5H which agrees with the above observed 
values of 8B. No decrease is observed at temperatures 
T > 3.4 K and near the beat nodes. 

This explanation is checked in a simple way. Usually 
the applied magnetic field H ramps much slower than the 
AC modulation field h, i.e. the full magnetic field H + h 



3 




3.6 4.0 4.4 3.58 3.60 



FIG. 2: (a) Pickup voltages divided by the modulation level 
for low and high modulation amplitude at 1.3 K. (b) DC 
susceptibility derived from magnetization measurements with 
Hall probes without field modulation. Graphs to the right 
show respective zooms. 
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FIG. 3: (a) Schematic representation of the hysteresis loop 
showing that the response to an applied AC modulation field 
is nonlinear, (b) Modulation level dependence of the normal- 
ized pickup voltage for characteristic magnetic fields at 1.3 K. 
(c) The same dependence at the antinode at T = 2.2 K and 
T = 4.2 K. Nonlinearities arise at a critical temperature of 
3.4 K. 



always oscillates around the quasi static offset field H 
and B{H + h) makes a loop in presence of hysteresis. 
If the ramp rates are changed in a way that dh/At and 
dH/At are the same, then the magnetic field sweeps only 
forward with small steps in the direction of the ramp. 
Under this condition B{H + h) never makes a loop as we 
go always along the hysteresis loop boundary and never 
inside. In this regime, which is usually not used, the lock- 
in amplifier does not measure the correct amplitude. We 
observed however in this regime the usual dHvA beat 
signal without notches. 
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FIG. 4: (a) Third harmoinc of the pickup voltage. The insert 
shows a zoom to the above discussed field range, (b) Out of 
phase part of the pickup voltage, (both at 2.5 G modulation 
level and 1.3 K). The insert shows a schematic representation 
of the hysteresis showing that the response to a sinusoidal 
field modulation becomes window shaped and is slightly phase 
shifted with respect to the input. 



The notches in the envelope of the first harmonic in- 
phase AC response are "compensated" by steeply rising 
higher harmonics at the same critical point. At the same 
field a phase shift appears in the pickup voltage. This 
means that the response to an AC modulation becomes 
extremely nonlinear in the presence of hysteresis. Fig. ^ 
shows the third harmonic (a) and the imaginary part (b) 
of the pickup voltage in a wide region of magnetic field at 
1.3 K. The amplitude is big in both curves only around 
the beat antinodes. In the regions of nodes the signal 
is about zero. The insert of fig. E^a) shows the 3'^'^ har- 
monic in the same field interval as in fig. [3 The compar- 
ison shows that the signal appears and disappears with 
a threshold character at the critical points of the transi- 
tion to the CDS. The insert in fig. EJb) shows schemat- 
ically that the response to an initially sinusoidal modu- 
lation field becomes highly distorted in the presence of 
a hysteresis loop. The response becomes more rectan- 
gular shaped like a window-function. This function can 
be composed, as it is well known, of odd harmonics of a 
sine. This explains that the third harmonic content of 
the pickup voltage is very high in the CDS. At magnetic 
fields without hysteresis the 3"^ harmonic is very small 
as the modulation level /i = 2.5 G is much less than the 
dHvA period of about 130 G at 3.6 T. 

The hysteresis shape and size is reconstructed in fig. |5l 
like e.g. in The response to the sinusoidal modula- 
tion of 2.5 G amplitude is calculated by adding all in and 
out of phase contributions up to the fifth harmonic. The 
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FIG. 5: Hysteresis shape reconstruction for modulation am- 
plitude 2.5 G for several characteristic magnetic fields at 
1.3 K. 



same procedure was applied at different positions along 
the dHvA oscillations. At the diamagnetic part (% < 0) 
of every dHvA period all harmonics vanish and only the 
in phase response persists. Hence, a line with negative 
slope is calculated. The same behavior is found for x > 
around a node with a line with positive slope. Whereas 
at the paramagnetic part in the region of the notches hys- 
teresis arises and its size is maximal at the dHvA beat 
antinode. 

The observed threshold behavior of the third harmonic 
and of the out of phase signal does not essentially change 
when the modulation level is varied providing h <^ AH. 
This behavior was observed in the frequency range from 
8 Hz up to about 200 Hz. Thus, these measurements 
offer a simple way to determine the Condon domain phase 
diagram. 

Moreover, hysteresis may explain the data measured by 
Plummer with an AC mutual inductance method at 
low modulation level. First Plummer thought that a new 
dHvA frequency was discovered. As the Fermi surface of 



beryllium was not consistent with this frequency, eddy 
currents at the rather high frequency of 100 Hz were in- 
voked 0,^^. A comparison of Plummer 's data with our 
results (see Fig. |2Ia), black curve) shows the similarity. 
The deep notches in the dHvA oscillation envelope seem 
indeed to be the result of a new frequency. We repeated 
the measurements at higher frequencies and found the 
same behavior. 

Hysteresis accompanies without doubt the appearance 
of the CDS resulting from the Pippard-Shoenberg insta- 
bility. For any sample shape, not only in plate-like sam- 
ples, in the CDS coexist two phases Bi and B2 with 
domain walls between them. The wall motion and pin- 
ning might, in principle, depend on the sample shape. 
Therefore, a plate-like sample was measured in the same 
T and H range. We found the same phase diagram for 
hysteresis formation in both samples. The hysteresis size, 
however, needs to be investigated more precisely as func- 
tion of temperature, magnetic field and for various n. 

Hysteresis itself is a result of interaction or pinning 
of domain walls with defects, impurities and the surface 
of the crystal. Moreover, the question of domain wall 
motion can not be neglected in many phenomena like e.g. 
acoustic wave propagation and absorption, and helicon 
waves Recently, the question of Condon domain 

wall motion was considered in detail theoretically [iSj. 
Unfortunately, an idealized model of a domain wall was 
used without taking into account the direct link with 
lattice deformation [14| and the metastable behavior due 
to hysteresis. 

In conclusion we have shown that hysteresis occurs in 
the dHvA effect under the conditions of the Condon do- 
main state. The observed hysteresis loop width is rather 
small, only a few Gauss. We have shown that the out of 
phase part and the third harmonic of the pickup voltage 
rise steeply when the magnetization becomes irreversible. 
This threshold character offers a simple and robust pos- 
sibility to measure a Condon domain phase diagram. 

We are grateful to I.D. Vaguer and I. Sheikin for fruit- 
ful discussions. 
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